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ABSTRACT 

The dynamics of detumbling a randomly spinning spacecraft 
using externally mounted, movable telescoping appendages are 
studied both analytically and numerically. Two types of telescoping 
appendages are considered: (a) where an end mass is mounted at the 

end of an (assumed) massless boom; and (b) where the appendage is 
assumed to consist of a uniformly distributed homogeneous mass throughout 
its length. From an application of Lyapunov's second method boom 
extension maneuvers can be determined to approach either of two desired 
final states: close to a zero inertial angular velocity state, and 

a final spin rate about only one of the principal axes. Recovery 
dynamics are evaluated analytically for the case of symmetrical deploy- 
ment. Numerical examination of other asymmetrical cases verifies the 
practicality of using movable appendages to recover a randomly tumbling 
spacecraft. 
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NOMENCLATURE 


hi » h 2 » h 3 


I P I 2 , I 3 

* * * 
Ip I 2 , I 3 

_* _* _* 
1 1 » 1 2 * I3 


(il ^ ^ (0 2 > U) ^ 


* 

4>o» 4)0 


« boom extension rate 

= components of the angular .momentum vector 
along the principal axes 

= constant value of h 3 for symmetrical 
deployment maneuver 

s instantaneous values of principal moments 
of inertia 

= hub principal moments of inertia 

= moments of inertia at the switching time T 3f 
in the recovery sequence to achieve final ' 
spin about the '3' axis 

. 3 

= 2 p c 

= time varying length of telescopic appendages 
= boom end mass 
= 2 me 2 
= ti me 

= kinetic energy 

= switching time in the recovery sequence to 
achieve final spin about the '3' axis 

= Lyapunov function 

= angular velocity components along the principal - 
axes 

= desir e --d final value of id 3 (u 3 ^) 

= mass density per unit boom length 

= phase angles appearing in the solutions for 
toi (t) , « 2 (t) and determined from conditions at 
t = 0, t = T 3 f , respectively 

= /b(t)dt 

= nutation angle 

= indicates time differentiation 

= indicates initial conditions 
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I. INTRODUCTION 

With the advent of prolonged manned missions in space and 

the possibility of in-orbit servicing and repairs to spacecraft 

there is an increased interest in studying methods that can be 

used to recover a spacecraft from an initial uncontrolled dynamic 

state. A recent paper examined methods of recovering spinning 

satellites to a flat-spin condition by using spin-up thrusters 

1 

and multiple combinations of thrusters. It was concluded that 

the use of such thrusters for the recovery operation are often 

limited by the weight and propellant capacity of the thruster 

system, and also the reliability problems associated with multiple 

thrusters in sequence. Kaplan has described an alternate recovery 

system which utilizes a movable-mass control device that is internal 

2 

to the spacecraft and can move along a fixed direction track. This 

device is activated upon initiation of tumble and is programmed via 

a control law to quickly stabilize motion about the major principal 

3 

axis. In a recent related paper it was concluded that the mass 
track should be placed as far as possible from the vehicle center of 
mass and be oriented parallel to the maximum inertia axis; in addition 
the performance of the control system can be improved through larger 
mass amplitudes along the track and also larger mass sizes. 

It is apparent that the location and displacement amplitude of 
any internal control mass will be limited by the physical dimensions 
of the space vehicle. 
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Externally movable appendages could allow for a greater range of 
location and displacement amplitudes of such a system; hov/ever, 
as the size of the appendages increases the flexibility problems 
associated with such structures would have to be considered. 

Of interest in this study will be the consideration of the 
detumbling dynamics of a spacecraft system with extensible boom- 
type appendages along the principal axes, The recovery maneuver from 
an initial tumble is designed to approach either of two final states: 
(!) close to a zero inertial angular velocity vector and (2) to approxi 
mate a final spin about a principal axis. (It is thought that small 
terminal residual angular rates could then be removed by temporarily 
activating on-board damping systems.) A key advantage of this type 
of system would be its potential reuse for subsequent detumbling 
recovery operations as the need arises. 

II. ANALYSIS 

A. General Considerations 

The dynamics of detumbling a randomly spinning spacecraft using 
externally mounted, movable telescoping appendages are studied both 
analytically and numerically. The appendages considered are of vary- 
ing length and could represent extensible booms or a tether connected 
to the main part of the spacecraft. Two types of telescoping append- 
ages are condidered: (a) the case where an end mass is mounted at 

the end of an assumed massless member (end mass moving) as shown in 
Fig.l (a);and (b) where the appendage is assumed to consist of a uni- 
formly distributed, homogeneous mass throughout its length (uniformly 
distributed mass moving) as shown in Fig. 1(b). 


The extensible boom type appendages are assumed to originate 
from the center of the hub along the three principal axes. The 
desired final states of the system considered are: (1) zero inertial 

angular velocity vector and (2) a final spin about one of the principal 
axes. The necessary conditions for(asymptotic)stability during the 
detumblirig sequences are determined using Lyaounov's second method. 

B. End Hass Moving 

1. Development of Kinetic Energy 

The configuration of the system, where the end masses are assumed 
to be attached to the end of massless rods along all three principal 
axes is shown in l : ig.l (a) .The end masses are assumed to be identical 
(i.e. m. ** m). The rotational kinetic energy of the system can be 
developed in terms of the hub inertias (I.j*) and boom lengths as: 

T « [{I* + 2mU 2 + & 3 )} Wj + {I 2 + 2m(& 3 + £i)}w 2 

2 2 2 2 , 2 *2 

+ { 1 3 + 2 m{&i + £ 2 )} w 3 * 2 m (& 1 + ft 2 * & 3 )] 

■ A 1 2 2 

Defining, 1 1 55 + 2m(& 2 + & 3 ) 

I 2 = I 2 + 

I 3 = Ig + 2m(£ 1 + it) , 

Eq. (1) can be rewritten as: 

1 ■ 2 2 ‘ 2 * 2 -2 , 2 

T = -g- [I o» 1 + I 2 <j 2 + I 3 oj 3 + 2m(£ 1 + i 2 + & 3 )] 


0 ) 

( 2 ) 

( 3 ) 


If the extension rates are assumed to be constant, Eq. (3) 
can be expressed: 


12 2 2 

T = g [I tij + I 2 w 2 + I 3 u 3 ] + non-negative const. (4) 


Here the moments of inertia are time varying as the length of the 

booms varies during extension. 

2. Achieve Zero Inertial Angular Rate 

(a) Lyapunov Function-Kinetic Energy 

The desired final state of the system is w. = 0. A suitable 

Lyapunov function, in the state variables io 2 and to 3 , is the 

system rotational kinetic energy which can be written as: 

2 2 2 

V = T = [I ^ + I w 2 + I ,w 3 ] + non-negative const. (5) 

The Lyapunov function, V, is positive definite in the state variables 
selected; for asymptotic stability V will now be examined. 

Differentiating Eq. (5) with respect to time, there results: 


V = \ (1^ + I 2 C0 2 + 1 3 lOg + 2I 1 (0 1 (0 1 + 2 1 2 + 2I 3 W 3 M 3 ) (6) 

Euler's equations of torque- free motion can be written in the follow- 
ing form: 

♦ ■ • 

hj s u 3 h 2 - w 2 h 3 = IjWj + I ico 1 (7a) 


h 2 = to 1 h 3 - ojghx = I 2 o > 2 + I 2 u> 2 

* ■ • • 
h 3 = M 2 hl - Mjh 2 = l 3 m 3 + 1^ 


(7b) 

(7c) 
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Multiplying Eq, (7a) by w 1( Eq. (7b) by w 2 , and Eq. (7c) by w 3 , 

and adding we obtain the following: 

. . , *2 '2 *2 

4- 12^20)2 + I 3 W 3 W 3 - - ( I jfDj + I 2 U 2 1 3^3 ) (8) 

Substituting Eq. ( 8 ) into Eq. ( 6 ), there results: 

v > - \ (ij«J + i 2 <4 + ijuj) 0) 

From Eq. (9), we conclude that V is a negative definite function in the 

state variables only if I 1 , I 2 , I 3 > 0 . 

Here it is seen that when the rotational kinetic energy is used 

as a Lyapunov function expressed in terms of small amplitudes of the 

inertial angular velocity components, that the necessary conditions 

for asymptotic stability are satisfied for positive constant boom 

extension rates and three orthogonally mounted sets of booms along 

the hub principal axes. This means that as time becomes extremely 

large (and boom lengths become infinite) it would be theoretically 

possible to achieve a zero inertial angular velocity state. (Of 

course, such a situation will, in practice, not occur due to finite 

length appendages and the presence of large amplitude rates for the 

general situation of an initial random tumble. However, it will be 

of interest to simulate how much of a random tumble could be removed 

by this process where stability is now considered in the global sense 

about w.. =0.) The selection of rotational kinetic energy as a Lyapunov 

3 

function has also been used by Edwards and Kaplan for the system treated 
in Ref. 2. 
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(b) Analytical Solution 

As a special case when the spin axis is an axis of symmetry 

(Ii n I 2 = I) during deployment, Eqs. (7) become: 

♦ 

hi + b(t) h 2 c 0 (10a) 

h 2 - b(t) hi = 0 (10b) 

h3 = h 0 = const (10c) 


where 

b (t) = Ia(t) - I(t) 


h o 



Introducing t = /b (t)dt Eqs. (10a) and (10b) reduce to, 


37- + h 2 = 0 (12a) 

- 111 = o (12b) 

which are in the standard form of the two dimensional harmonic 
oscillator. The solutions to Eqs. (12) can be written as 

h i ( t) = q* cos t ~ q q cos (/J b(t)dt + ij' 0 ) (13) 

h 2 (t) = q* sin t qj sin (/* b(t)dt + *J) (14) 

The solutions given by Eqs. (13) and (14) are’ identical with those 
previously given by Hughes as a special case of his approximate 
analytical solution for the motion during deployment of a spacecraft 

5 

with telescoping booms where the nutation angle remains small. 


Wo apply this solution to the symmetrical configuration of 
Fig. 1(a) where the moments of inertia about the principal axes 
are expressed, for the case of a uniform extension rate, c, along 
all three axes, by; 


I x = I -Mm % - I + 2Pt 

* g * 2 

I„ = I + 4m* = I + 2Pt 


1 3 * I g + W* = 1 3 + 2Pt Z 


and P a 2mc 


Using Eqs. (15) in Eq, (11) we obtain; 


- ^r+-p 
2 


where = 71*72 P' and d 2 = /I - */ 2P' (17) 

Introducing Eq. (16) in Eqs. (13) and (14), and after performing the 
integration, the solutions for the transverse angular velocities are 
obtained as: 


» x (t) = 


* n 0 i -l t 

q 0 cos[ w { ^tan ^ 

I* + 2 Pt* 


i -if *„ 

HT ta " g- )+ V'o] 


% C 5F < 3T ‘a"'’ 3T ' k ta " 

I* + 2Pt^ 


1 t , *n 

HT } + ^ 
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where q 0 and <ij are determined from the initial conditions. 

From Eq. (10c), the angular rate about the *3' axis is written, 

I* »,(0) 

u 3 (t) = rf"“2PF (20) 

we observe here for large values of t, the solutions for the angular 
velocities lead to the form: 

w 4 (t) = const/(l* + 2Pt 2 ), i = 1,2,3 (21) 

This equation indicates that the magnitudes of the angular velocities 
decrease during extension of the appendages, with the square of the 
elapsed time. 

3. Achieve Final Spin About One of the Principal Axes 
(a) Lyapunov Function-Modified Kinetic Energy 
The desired final state of the system is: uj = 0, io 2 = 0 and 
w 3 = u 3.f = Using the state variables: <o 2 » and <03 - n» the 

Lyapunov function is defined as the modified rotational kinetic 
energy, which can be written as: 

V - \ [Ijuj? + hi + h (“3 - (22) 

Here V is positive definite in the state variables selected. Differ- 
entiating Eq. (22) with respect to time, there results: 

V s 2 [i^l + V J 2 * “ n ) 2 

+ + 2 + 2 1 (tu 3 - fl ) £03 ] 


(23) 


Using Eq. (3) in Eq. (23), wg obtain, 


• i*2 *2 • 2 I 2* • 

V a ~ + l 2 w?. + l 3 w 3 )+gn I 3 - (I3W3 + 

For symmetry about tho *3' axis during extension: 



( 24 ) 


h 3 = I 3 W 3 + I 3 U> 3 ” 0 

Eq. (25) is used in Eq. (24) to obtain: 

. 1 . 2 • 2 * 2 2 

V « - 7 [IjWj + 1 2^2 + I 3 ( w 3 " t? )] 

After rewriting Eq. (25) in terms of the state variables, 

V = - •]? [i^i + i 2 u >2 + l 3 (w 3 -n) 3 - i 3 n (n > 3 - n) 


(25) 


(26) 


(27) 


Also from Eq. (25), the solution for w 3 (t) is given by, 

w 3 (t) = Igu 3 (0)/I 3 (28) 

We conclude from Eq. (27) that V is negative definite in the state 
variables only if: 

» * # 

w 3 i fi, Ip I 2 > 0, and I 3 > 0 for w 3 > n 

Thus for the case where a spin about one of the principal axes 
is a desired final condition, a modified form of the kinetic energy 
can be used as a Lyapunov function. 
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Here the final state can be achieved by extending all telescoping 
booir.s until the desired spin rate is reached and then continuing 
the extension of the set of booms along the nominal spin axis until 
the transverse components of angular velocity reach an acceptably 
small amplitude (within the limitations of boom length), It should 
be noted that if we allow w 3 < n and I 3 f 0> there will be a 
difference in sign between the third and fourth terms in Eq. (27). 
(b) Analytical Solution 

The time at which w 3 « to 3 ^ = n will be denoted by T 3 ^ > 

At t = T 3f> 

■If ^ 

I, » I 2 ■ I* + 2P (T 3f ) 2 (29) 

-* 2 

I, - lj + 2P (T Jf ) (30) 


For t < T 3f the solutions for the angular velocities can he obtained 
from Eqs. (18), (19) and (20). 

For t > T, 


3f > 


.* 

■ 

•3 

.* 


I 3 = const = I* + 2P (T 3 )' 


where 


•if 2 2 2 

i « r + p(t 3f ) + pt = i f + pt 


If - I + P(T 3f ) 


From Eq. (11), and using Eqs. (31) and (33), we obtain 


_* 


b(t) ■ % 1 - 11 
f 


(31) 

(32) 

(33) 


(34) 
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Introducing Eq. (34) into Eqs. (13) and (14), the solutions for 
the angular velocities for t > T 3 are, 

, r*~ U - ' 

q 0 c °s Cw 3 t^i* p {tan )-tan (-p= 

' f f 


~) t+T 3 ]+t 0 ] 
I*/p f 


to x (t) = 


(I* + Pt 2 ) 


(35) 




% sinC % [ 7=T {ta,, ' I( '^r )J-t+T, f >* 0 ] 


to 2 (t) * 


(I* + Pt2) 


(36) 


and from Eq. (10c), 


w 3 (t) « w a , ~ const (37) 

"f 

Here q Q and ^ are to be determined from Eqs. (18) and (19) at 
t = T 3 ^ and should not be confused with q$ and i|j* which are determined 
at t * 0. 

For large values of t, Eqs. (35) and (36) reduce to the form, 


q Q cos (w 3 ^ x const x t + const) 


“l = 


I* + Pt? 


(02 = 


q 0 sin (a> 3 ^ x const x t + const) 
“ I* + .Pt* 

f 


(38) 


(39) 
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The above two equations indicate that the frequency of oscillation 
approaches a constant value and the magnitude of the oscillation 
decreases with the square of the elapsed time. 

The time, t = T 3f , at which the extension of the booms along the 
1 1 ' and '2' axes are stopped can be determined from « 0, yielding 
the result: 



1 

2c 



/ m 3 (°> 

l ) 


(40) 


C. Uniformly Distributed Mass Moving 
1. Achieve Zero Inertial Angular Rate 
(a) Analytical Solution 

The desired final state of the system is w. = 0. The booms con- 
sidered are assumed to have a uniformly distributed mass ( p) along 
their lengths. The same procedure as adopted in the case of the 
moving end masses can be applied here to obtain the solutions for 
the angular velocities. Here we present only the final results. 

The solutions for the angular velocities are given by: 


(Oi (t) 


w 2 (t) 


w 3 (t) 


q* cos { ^ b(t) dt + } 

I* + § Kt 3 

q* sin { b(t) dt + \J>* } 

o ij o 

r*+ I Kt 3 

I 3 w 3 (0) 

"* 2 3 

I 3 + 3 Kt 


and K = 2 pc 3 


(41) 


(42) 
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In Eqs. (41) and (42), 


/t b(t) dt . Zl r _i_ an{ , (d 8 ± t ,) 2 >}+ ill d .i } 

ft tN L CJ ? J9 i.j.9 ■> T j? /"TT" 1 ' 11 x ^ yg — ■> 

3 


6d| 


d| - d 3 t+t 2 


df^“ 


6d 


1 r „ jdSL±tl j „ 


4 u 4 


-T 2 ~: tan 1 { • ~r= } ] 

vr d h rr 


(44) 


3 31 3 

^3 a snd dt^ 


31, 

2K 


(45) 


2. Achieve Final Spin About One of the Principal Axes 
(a) Analytical Solution 

The desired final state of the system is ° 0, to 2 = 0 and 
w 3 = w 3 ^ = n. For t i T 3 ^ s the solutions for the angular velocities 
can be obtained from Eqs. (41), (42) and (43). For t > T 3 ^, the 
solutions for the angular velocities can be obtained as: 


Mjtt) 


w 2 (t) 

to 3 (t) 

£ 

where I, 
f 


q cos{ b(t) dt + f } 


I* + ^ K t 3 

q s i n { b(t) dt + ¥ } 


* 1 3 

I + 4r K t 
f J 


to 9 = const 
3 f 

,* + | (T, f >' 


(46) 


(47) 


(48) 

(49) 
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In Eqs , (46) and (47), 


3 I 

/ t b(t)dt = w 3 [ — 

T 3 f " 

f 


1 (d s + t) 2 

6d2 * n { c! 2 - d t + t 2 } 
5 5 5 




tan 


-l 


2t - d _ 

5 

d. rT~ 

D 


} ) -t + T, 3 
3 f 


(50) 


_* 

I. 



K (T 3 ) 3 
3 f 



(51) 

(52) 


The time T , at which the booms along the and '2' axes are stopped, 
can be obtained as: 


3 I3 

T = [ ( 

3 L 2K v 


w 


,(0) -w 3j 


3* ' ” a f % n i/a 


) 3 


(53) 


III. NUMERICAL RESULTS 


A. End Mass Moving 

A typical detumbling maneuver for an initially uncontrolled 
spacecraft is illustrated in Figs. 2 and 3. In this example because 
of symmetry the uncontrolled torque- free nutation (Fig. 2) can be 
theoretically predicted. This motion is also represented by the 
dotted curves corresponding to a zero boom extension rate (c. a 0) 
in Figs. 3. The effect of extension rate on the recovery is illustrated 
in Figs. 3(a) - (c). For small extension rates (up to 1 ft/sec) the 
oscillatory nature of the transverse motion is not removed until after 
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the first cycle; the advantage of considering higher extension rates 
(at the expense of on-board power) for an initial fast tumbling is 
apparent. It should be noted that at a given time in these figures 
different boom lengths are represented according to the extension 
rate. For example, with an extension rate of 4 ft/sec after 60 ft, 
of extension along all three principal axes the angular velocity 
components have been reduced by more than a factor of 10 and, if 
240 ft. of boom could be extended, by a factor of over 300. Removal 
of this residual angular velocity could then be achieved by temporarily 
activating on-board damping devices. Then the appendages could be 
retracted and would be ready for subsequent reuse as necessary. 

Numerical examination of other cases for asymmetrical hubs 
also verifies the practicality of using movable appendages for the 
initial detumbling of randomly spinning spacecraft (Figs. 4(a) and 
(b ) ) . The numerical simulation results for an asymmetrical spacecraft 
are compared with the closed form solution for a symmetrical extension 
and it is observed that the closed form solutions are only applicable 
when the asymmetry is small. Although the moving end mass system is 
effective in reducing the magnitude of all components of the angular 
velocities, for the masses considered here it will not effectively 
reduce the initial nutation angle during the response time simulated. 

For the symmetrical deployment shown in Fig. 4 the initial nutation 

y 

angle of 27.89 degrees is ihaintained, whereas for the asymmetrical 

case, an initial nutation angle of 29.09 degrees is increased within 

3.5 seconds to 36.95 degrees and maintained up to the 60 seconds simulated. 
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The nutation angle is defined here as the angle between the nominal 
symmetry ('3') axis and the total angular momentum vector. 

Fig. 5 illustrates a recovery maneuver which would result 

in a final spin about the '3' body axis with a small transverse 

residual. The booms are extended so that the modified rotational 

kinetic energy is positive definite and its total time derivative 

is negative definite during the maneuver. All booms are extended 

until T 3 ^ at which time io 3 = u) 3 ^. Then, only booms along the + '3' 

axis are extended to reduce the transverse residual components. At 

the switching time, T a ., the effect of any lag in the system has not 

f 

been considered here but should be considered in the design and per- 
formance of the actual system. From the figure it can be seen that 
the frequency of the response of the transverse components of the 
angular velocity is essentially constant in the terminal part of the 
maneuver (consistent with the discussion in connection with Eqs. (38) 
and (39)). 

A comparison of the recovery maneuver of an asymmetrical spacecraft 
with that of a symmetrical spacecraft to achieve a final spin along the 
'3' axis is shown in Figs. 6(a) and (b). The calculated T 3f for the 
symmetrical spacecraft is used for stopping the booms along the 'V 
and '2' axes. It is observed that using this logic the final w 3 

f 

reaches a lower value (1.8 rad/sec) when compared with the desired 
final value (2.0 rad/sec). Also we notice from Fig. 6(a), the response 
of iOj(t) for the asymmetrical case differs from that of the symmetrical 
case. 
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This is due to the increase in the. order of the system equations 

for the asymmetrical extension (i.e. - three first order differential 

equations must now be considered). It should be pointed out that 

after T 3 ^., for the asymmetrical case, the time response of w 3 is not 

exactly a straight line as apparently indicated in Fig. 6(a) but also 

consists of small amplitude oscillations superimposed about this 

straight line solution. For larger asymmetries this oscillation 

would become apparent within the plotting scale shown and the 

difference between w, achieved and desired would also increase 

3 f 

using the open loop control logic of switching the extension sequence 

at a pre-set T 3 . It should also be noted that the nutation angle, 
f 

for both cases, does not vary greatly after T 3 ^, 

B. Uniformly Distributed Mass Moving 
For the case of a spacecraft with a uniformly distributed 
mass along the boom lengths a typical detumbling maneuver is illustrated 
in Figs. 7(a) and (b) for a symmetrical and asymmetrical hub, respectively. 
If 240 ft. of boom could be extended the total mass of each appendage 
would be about 1 slug. The effect of the increased appendage mass 
(when compared with the case of Figs. 3 and 4) i s immediately apparent 
by noting the extremely small amplitude of the residual angular velocity - 
of the order 10“^ rad/sec. The nutation angle behavior for both cases 
of deployment is essentially the same as already described in connection 
with Fig. 4. 

A representative recovery operation to achieve a final spin about 
the '3' axis with this system is shown in Figs. 8(a) and (b). 
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A comparison with the results of Figs. 6(a) and (b) clearly demonstrates 
the effectiveness of the heavier boom after the switching time. In 
addition, after approximately 40 seconds a drastic reduction in the 
nutation angle is also apparent - a result not achieved by the lighter 
end mass system. Despite the obvious advantages of using a more 
massive appendage there are, however, two practical considerations 
that have been ignored: (1) the increased weight of the payload 

package and (2) the effect of flexibility on the assumed rigid body 
dynamics simulated here. 

IV. CONCLUDING COMMENTS 

As an application for spacecraft rescue and recovery, when booms 
are extended along all the principal axes to detumble a symmetrical 
spacecraft, exact closed form analytical solutions can be obtained 
for all three angular velocities of the spacecraft. 

Boom extension maneuvers can be determined to approach either 
of the two desired final states using Lyapunov's second method. The 
conclusions are that: (1) as time becomes extremely large (and boom 

lengths become infinite) it would be theoretically possible to achieve 
a zero inertial angular velocity state (of course, such a situation will, 
in practice, not occur due to finite length booms); (2) the final spin 
about one of the principal axes can be achieved by extending all tele- 
scoping booms until the dpsired spin rate is reached and then continuing 
the extension of the set of booms along the nominal spin axis until the 
transverse components of angular velocity reach an acceptably small 
amplitude. 
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Numerical examination of other cases for asymmetrical hubs 
also verifies the practicality of using movable appendages for 
the initial detumbling of randomly spinning spacecraft, 

An advantage of the telescoping system as used in the recovery 
of tumbling spacecraft is its potential reuse. The booms can be 
retracted at the end of each recovery operation once the snail 
residual angular velocity components have been removed by temporarily 
activating on-board damping devices. The constraints on such a system 
are: (!) the limitations on the extension rate, and boom lengths that 
are practicable; (2) the limitations on the rate of initial tumble 
that could be handled by the system without compromising its structural 
integrity, 
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FIG. 2. DYNAMICS OF UNCONTROLLED MOTION 
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